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TECHNICAL MEMORANDUM 1441 


THE EFFECT OF FREE- STREAM TURBULENCE ON 

HEAT TRANSFER FROM A FLAT PLATE* 

By Sugao Sugawara, Takashi Sato, 
Hiroyasu Komatsu, and Hiroichi Osaka 


The effect of initial turbulence on the surface heat transfer was 
investigated in the case of a smooth flat plate. The initial turbulence 

was introduced by a screen, and the percentage turbulence was 

measured accurately by a hot-wire anemometer in the boundary layer and 
in the main flow. On the other hand, the change of the temperature dif- 
ference between air and a plate, during the plate's being cooled, was 
measured in the same condition, and from its result the local heat- 
transfer coefficient on the surface was calculated. 

Thus, we were able to make clear the relation between initial tur- 
bulence and heat transfer. Surface heat transfer becomes better with 
the enlargement of the initial turbulence, and the increasing of heat- 
transfer coefficient was remarkable within the range of small initial 
turbulence. In the rage of larger initial turbulence than 7 to 8 per- 
cent, the local heat- transfer coefficient increased no more, and its 
value was 55 percent larger compared with the case of the smallest ini- 
tial turbulence. 


1. INTRODUCTION 

Surface heat transfer Is a phenomenon of heat exchange through the 
boundary layer surrounding the surface of a solid body, it is natural, 
therefore, that the factors affecting the condition of boundary layer 
also Influence heat- transfer phenomena. It is considered that turbulence 
existing In the main flow will Influence the boundary layer and, there- 
fore, surface heat transfer. Many investigations of surface heat trans- 
fer have been conducted in which very little or no turbulence exists in 
the main flow, or without consideration of any turbulence that may be 
present. However, in many practical situations, comparatively strong 
turbulence is present and, therefore, should be taken into account. 


* "Shuryubu midare ga heiban nets udentatsu ni oyobosu likyo ni 
tsuite," Jour. Jap. Soc. Mech. Eng., vol. 19, no. 18, 1953, pp. 18-25. 
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Experimental results that were made to clarify the effect of such 
turbulence on heat transfer at the surface of a smooth flat plate are 
reported herein. The turbulence was generated by screens in the wind 

tunnel, and the turbulence percentage -y^T^/u was measured accurately 
by a hot-wire anemometer both in the boundary layer and the free stre am . 


2. EXPERIMENTAL APPARATUS AND METHOD 


i. 


The main apparatus as shown in figure 2 was the same as that de- 
scribed in reference 1. An apparatus to generate various degrees of 
turbulence was added as shown in figure 1. It consists of a screen S 
behind the tunnel inlet T , built by connecting a number of circu lar 
rods in a square array. The important dimensions of the rods are given 
in the following table: 


TABLE 1. - SCREEN DIMENSIONS, mm 


Screen 

A 

B 

C 

Rod diameter 

6 

3 

1 

Rod pitch 

28 

13 

7 


The turbulence generated was the strongest when using arrangement A 
and the weakest with C . A smooth flat plate P was suspended in a 
duct M. Since the turbulence induced by the screen decreases with in- 
creasing distance from the screen, the turbulence near the measuring 
plate P can be adjusted by varying the relative position of the screen 
and the plate. This screen- to-plate distance is changed by an insert 
duct D between the Bcreen and the plate which had a length of either 
50 or 10 centimeters. Experiments were carried out by combining two 
lengths of D, or without D, with three kinds of S. 


ii. [sic] Heat Transfer Coefficient Measurements 

The method of measurement was the same as described in reference 1. 
The local heat- transfer coefficient a is calculated from the relation 

a = Bey s/2 _ (l) 

where c, y, and s are specific heat, specific weight, and thickness 
of a plate, respectively. The constant B multiplies the time t 
in the equation of temperature difference between the plate and air 
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0 , as defined by 


0 = Ae‘ Bt (2) 

B was determined from the relation between 0 and t which was meas- 
ured directly by two millivolt meters. Local heat transfer coefficients 
were measured at twelve locations along the plate, and two values of free- 
stream velocity were used (14 m/sec and 7 m/sec). 


iii. Measurement of Turbulence Strength 

(a) Frequency characteristics of amplifier used to measure turbu- 
lence . - A hot-wire anemometer was used for the measurement of turbu- 
lence. With this type of sensing element, the frequency response is 
limited by the heat capacity of the hot wire. That is, the temperature 
variations of the hot wire which are read as voltage variations between 
the terminals do not follow the velocity variations of the flowing fluid, 
and both amplitude attenuation and phase shift are present. Moreover, 
both of these effects become greater as the frequency increa ses . Fro m 

reference 2, the amplitude decrease is proportional to l / sjl + (Ms) 2 , 
and the phase shift is proportional_to tan _ 1 Mu, where o> = 2 *f, with f 
as the frequency and M = 4.2 pA 2 c(T - T())/i 2 r proportional to the heat 
capacity of the hot wire, p, A, c, and r are specific resistance, 
cross-sectional area, specific heat, and resistance of the hot wire, 
respectively, T and Tq are mean temperatures of the hot wire and 
fluid, respectively, and i is the current through the hot wire. 

Such amplitude attenuations and phase shifts must be compensated 
since turbulence is distributed over a frequency spectrum. In order to 
accomplish this compensation, the characteristics of the amplifier must 
be the reverse of those of the hot wire. 

The compensating method used in the amplifier utilized negative 
feedback as described in reference 3 . The circuit diagram is shown in 
figure 2. Negative feedback is accomplished from the plate of V 3 to 
the cathode of V 2 > 

G(Ei + EqP) = Eo 

where G is the _ amplification of V 2 and V 3 , 0 is the amo unt of nega- 
tive feed-back, Erj_ is the input voltage of Vg, and Eq is output volt- 
ate of V e . Therefore, 


* 


®o/ B i = A = 


1 - PG 
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where A Is the actual amplification including negative feedback. 

Specially., when ( G^{ > > 1, 

A s - i (3) 

P 

Referring to figure 2, the bombined impedance z of Rj and Cj is 

• R 1 . rr 

z 1 + joiCiRi J . v- 1 

and if we choose Ro»|z| , 

P = z/(R 0 + z) s? z/R 0 

Therefore, from equation (3), 

A = - B 0 /z = - R 0 (l + JCDCiRiJ/Rq (4) 

If we choose the values Cp and Rp as CpRp = M, we can compensate for 
the decrease of amplitude and phase shift of the hot wire. It will be 
noted that the characteristics of the amplifier must be as flat as pos- 
sible when negative feedback is not applied. “ 

A vacuum thermocouple was used for the measurement of the alterna- 
tive component (turbulent component) of the output voltage. C is a 
blocking condensor for the constant (nonvariable component) voltage. 

Voltage induced in the T-circuit is proportional jto -yu 2 , where 
u is variation of velocity (turbulence); a bar over a symbol means an 
average with respect to time. 

Since the diameter of the hot wire used in our experiment was 30 
microns and its length was about 10 millimeters, the current passing 
through a hot wire was adjusted to 300 milliamperes. (The average tem- 
perature of the hot wire was about 170° C for a mean flow velocity of 14 
m/sec, so that M= 0.023 sec.) 

Figure 3 shows the frequency characteristics of our system where 
the broken line is the theoretical result of-combining a hot wire with 
an amplifier having a perfectly flat characteristic; the circles show 
the results of combining a hot wire with an amplifier without the nega- 
tive feedback circuit. The double circles show the output of the T- 
circuit when a hot-wire anemometer was connected with our amplifier using 
negative feedback. As shown in figure 3, the curve for" the double cir- » 

cles is satisfactorily flat between 20 and 7000 cycles. Figure 4 shows 
the calibration equipment for determining frequency characteristics. 
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The calibration for the absolute value of -/u 2 was carried out by- 
using a dynamic speaker. A hot-wire probe was fixed to the moving part 
of the speaker and the amplitude was measured by a traversing micrometer. 
The frequency was checked by a tuning fork and standard frequency meter. 

(b) Effect of mean velocity . - The sensitivity of a hot-wire may be 
affected by the change of mean velocity. The relation between the voltage 
E of the direct current that appears between the hot-wire terminals, and 
the mean velocity U is not in a straight-line relation as shown in fig- 
ure 5. This means that the sensitivity of a hot wire is a function of 
the mean velocity. Namely, its sensitivity must be proportional to 
dE/dU at the point of U]_ when the mean velocity is U]_. Therefore, we 
can make a correction by multiplying the ratio dU/dE by the value at 
some standard condition where the calibration was made. The relation 
between U and E depends also on the fluid temperature as shown in 
figure 5. The broken line in figure 5 shows the value of dU/dE at 18° C. 
The correction for temperature is necessary, of course, but the change of 
dE/dU is not as large as the change of temperature, since the curves for 
variable temperature are approximately parallel. 


Figure 6 shows the circuit diagram for the apparatus used to measure 
the mean velocity. In principle it is a bridge circuit that utilizes a 
difference method. R is a hot wire and the voltage difference appears 
between the terminals of R and R]_. (Rj_ is a constant resistance that 
is scarcely affected by temperature changes.) The currents of the two 
circuits are kept constant by adjusting Rq and Rg. 


Utilizing the above mentioned, we can measure the value of 
which is the percentage of turbulence . 



3. DISCUSSION OF RESULTS 

i. Turbulence and Velocity Distributions in the Boundary Layer 

Figures 7 and 8 show two examples of the experimental results. The 
distance from the plate surface is y, and f, m, and r indicate plate 
position as shown in figure 1. (Front edge of a plate, f; middle of a 
plate, m; near rear edge of a plate, r.) Figure 7 illustrates the results 
when no screens are used but with the auxiliary duct D. Consequently, 
the turbulence is comparatively small. Figure 8 shows the result of using 
screen A but without duct D, and is indicative of the effect of strong 
turbulence in the main flow. Broken lines in both figures show ^a/u 2 /U, 
with U the local mean velocity at y. The full lines show -yu 2 /Uo 
where Uq is the constant mean velocity of the main flow (approximately 
14 m/sec in the case of fig. 7 and 12 m/sec in the case of fig. 8). At 
the positions m and r, the boundary layers were turbulent in both 
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cases. The value of -^/G2/u reaches a maximum value of one-third or 
one-fourth of the boundary- layer thickness, and its value is larger in 
the case of figure 8 than in figure 7. The thickness of the boundary 
layer is considerably greater in the case of figure 8. By comparing 
those two figures, we can recognize that the percent of turbulence in the 
boundary layer becomes larger through the increase of the turbulence 
level of the main flow; at the same time, the thickness of the l ami nar 
sublayer becomes smaller. (This is also apparent in figs . 9 and 10.) 
Figures 9 and 10 show the time - mean velocity distribution in the bound- 
ary layer in the same cases as figures 7 and 8. Open circles show the 
results for the case of figure 7, and solid circles show those of figure 
8. The thickness of the boundary layer 8 was determined as the value 
of y, where the velocity was 99.5 percent of the main-flow velocity. 


ii. Heat-Transfer Results 

Figures 11 and 12 show heat- transfer results. Those figures show 
the relation between Nusselt number Hu and Reynolds number Re in a 
few cases of variable turbulence, where - — 

Nu = ax/X Re = Uqx/v 

and a is the local heat- transfer coefficient at a distance x from the 
front edge. 

Figure 11 shows the case for screen A with a duct length D of 
50 centimeters; the influence of the mean velocity of the main flow is 
shown. There is no appreciable effect of mean velocity on the relation 
between Re and Hu, which is similar to the case of low main-stream 
turbulence. Laminar heat transfer is shown as a square-root relation, 

Hu = 0.425 Re 0,5 . In the laminar boundary- layer range the effect of tur- 
bulence is not so remarkable, but the transition Re from laminar to 
turbulent flow decreases as the turbulence in the main flow increases. 

Over the entire range of main- flow turbulence studied,_ the relation be- 
tween Hu and Re can be expressed as follows for the range of turbulent 
heat transfer: 

Hu = KRe 0 - 8 (5) 

This slope is the same as the case that includes very small turbulence in 
main flow, i.e., the straight lines that represent the relation of tur- 
bulent heat transfer are parallel to each other as shown in figures 11 and 
12. The value of Hu increases according to the increase of main-stream 
turbulence. In other words, the value of K in equation (5) increases 
with the increase of turbulence in the main flow. Figure 13 shows the 
relation between percent of main- stream turbulence and the increase of 
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Nu. The abscissa is (Nu - Nuq)/Nuq, where Nuq is the value of Nu in 
the case of the smallest turbulence (percentage turbulence was 0.37 per- 
cent) in our experiments and the relation between Nuq and Re is shown 
as follows : 

Nu 0 = 0.0194 Re 0 - 8 (6) 

The coordinate of figure 13 is the percentage turbulence in the main flow 
near the front edge of the plate (point f in fig. l) . The main flow 
turbulence decreases in the flow direction in this case, so it is diffi- 
cult to determine the position of turbulence to be used for such figures. 
It Is noted that the turbulence near the transition point affects mainly 
the turbulent boundary layer heat transfer. It is difficult, however, to 
determine the transition point exactly and moreover, the transition point 
existed comparatively near the front edge in these experiments since the 
form of the front edge was flat and not sharp. Therefore, the turbulence 
near the front edge Is used in figure 13 for convenience. The following 
table shows the result of decreasing the turbulence In the main flow. 


TABLE 2 


Length of duct D, cm 

None 

None 

50 

None 

Kind of screen 

None 

A 

A 

B 

Decrease of turbulence in main 

flow. 

percent 

Position f 

TnACA note: 

0.75 

7.42 

2.05 

3.54 

Position m 

see fig. l] 

0.72 

4.85 

1.53 

2.68 

Position r 


0.75 

3.58 

1.32 

1.78 


As shown In figure 13, Nu increases with increasing turbulence, and the 
rate of increase is larger in the low turbulence range than in the high 
turbulence range. In the range larger than 7 to 8 percent turbulence, 

Nu becomes almost constant. The maximum increase of Nu was 55 percent. 

The inclination of such changes is very similar to the results re- 
ported by Comings, et al. (ref. 4) for the case of flow perpendicular to 
a circular cylinder. In their experimental results, Nu did not change 
any more in the range between 7 and 22 percent. We cannot, however, 
quantitatively compare our results with theirs, since not only are the 
geometries different, but also the then- reported mean heat- transfer coeffi- 
cients, including the laminar heat- transfer range. 
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4. CONCLUSIONS 

In order to study the effect of main-stream turbulence on the sur- 
face heat transfer in the case of forced convection along the flat plate, 
experiments were carried out by changing the main-stream turbulence 
through the use of screens and settling sections. The results obtained 
are summarized as follows: 

i . The transition Reynolds number from laminar to turbulent flow- 
decreases as the main-stream turbulence level increases. 

ii . Increasing the turbulence of the main flow, increases the tur- 
bulence level in the turbulent boundary layer, and decreases the thick- 
ness of the sublayer. This causes an increase in the heat- transfer 
coefficient. 

iii. In the range of laminar heat transfer, the effect of turbulence 
in main flow was not great. In the range of turbulent heat transfer, the 
heat- transfer coefficient increases according to the increase of tur- 
bulence, but Nu was still proportional to Re^-S. 

iv. Nu increases sharply in the smaller turbulence range, and the 
increase becomes gradual in the range of larger turbulence. Nu is 
almost constant in the range of turbulence larger than 7 to 8 percent; - 
The maximum value of Nu was 55 percent larger than the minimum value 
of Nu. 
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Figure 2. - Circuit diagram of amplifier. 20 mV 
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Figure 3. - Frequency characteristics. 
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Figure 7> - Turbulence distribution in boundary layer. (No screen, no 
auxiliary duct.) 


Figure 8. - Turbulence distribution in the boundary layer. (Screen A 
no auxiliary duct.) 






Figure 9. - Velocity profile in boundary layer at the middle of plate. 
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Figure 11. - Eolation "between Nu/Be and Re (screen A; duct, 50 era.) 














